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We investigate chiral symmetry breaking in strong magnetic helds at finite temperature and 
densities in a 3 flavor Nambu Jona Lasinio (NJL) model including the Kobayashi Maskawa t- 
Hooft (KMT) determinant term, using an explicit structure for the ground state in terms of quark 
antiquark condensates. The mass gap equations are solved self consistently and are used to compute 
f~>) ■ the thermodynamic potential. We also derive the equation of state for strange quark matter in the 

£f>^ ' presence of strong magnetic fields which could be relevant for proto-neutron stars. 
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I. INTRODUCTION 

The structure of QCD vacuum and its modification under extreme environment has been a major theoretical and 
experimental challenge in current physics (lj. In particular, it is interesting to study the modification of the structure 
of ground state at high temperature and/or high baryon densities as related to the nonperturbative aspects of QCD. 
This is important not only from a theoretical point of view, but also for many applications to problems of quark 
gluon plasma that could be copiously produced in relativistic heavy ion collisions as well as for the ultra dense cold 
nuclear/quark matter which could be present in the interior of compact stellar objects like neutron stars. In addition to 
hot and dense QCD, the effect of strong magnetic field on QCD vacuum structure has attracted recent attention. This 
is motivated by the possibility of creating ultra strong magnetic fields in non central collisions at RHIC and LHC. 
The strengths of the magnetic fields are estimated to be of hadronic scale @, Q of the order of eB ~ 2mJ (toJ ~ 10 18 
Gauss) at RHIC, to about eB ~ 15m£ at LHC @|. 

There have been recent calculations both analytic as well as with lattice simulations, which indicate that QCD 
phase digram is affected by strong magnetic fields • One of the interesting findings has been the chiral magnetic 
effect. Here an electric current of quarks along the magnetic field axis is generated if the densities of left and right 
handed quarks are not equal. At high temperatures and in presence of magnetic field such a current can be produced 
locally. The phase structure of dense matter in presence of magnetic field along with a non zero chiral density has 
recently been investigated for two flavor PNJL model for high temperatures relevant for RHIC and LHC 7]. There 
have also been many investigations to look into the vacuum structure of QCD and it has been recognized that the 
strong magnetic field acts as a catalyser of chiral symmetry breaking [HHI . 

In the context of cold dense matter, compact stars can be strongly magnetized. Neutron star observations indicate 
the magnetic field to be of the order of 10 12 -10 13 Gauss at the surface of ordinary pulsars [l2| . Further, the magnetars 
which are strongly magnetized neutron stars, may have even stronger magnetic fields of the order of 10 15 — 10 16 Gauss 
[131 - 4191 ] - Physical upper limit on the magnetic field in a gravitationally bound star is 10 18 Gauss that is obtained by 
comparing the magnetic and gravitational energies using virial theorem [l2[ ■ This limit could be higher for self bound 
objects like quark stars [20]. Since the magnetic field strengths are of the order of QCD scale, this can affect both 
the thermodynamic as well as hydrodynamics of such magnetized matter [2l| ■ The effects of magnetic field on the 
equation of state have been recently studied in Nambu Jona Lasinio model at zero temperature for three flavors and 
the equation of state has been computed for the cold quark matter [22| . It will be interesting to consider the finite 
temperature effects on such equation of state, which could be of relevance for proto-neutron stars. 

We had earlier considered a variational approach to study chiral symmetry breaking as well as color superconduc- 
tivity in hot and dense matter with an explicit structure for the 'ground state' (23l - [25| . The calculations were carried 
out within NJL model with minimization of free energy density to decide which condensate will exist at what density 
and/or temperature. A nice feature of the approach is that the four component quark field operator in the chiral 
symmetry broken phase gets determined from the vacuum structure. In the present work, we aim to investigate how 
the vacuum structure in the context of chiral symmetry breaking gets modified in the presence of magnetic field. 

We organize the paper as follows. In section II, we discuss the ansatz state with quark antiquark pairs in the 
presence of a magnetic field. We then generalize such a state to include the effects of temperature and density. In 
section III, we consider the 3 flavor NJL model along with the so called KMT term-the six fermion determinant 
interaction term which breaks U(l) axial symmetry as in QCD. We use this Hamiltonian and calculate its expectation 
value with respect to the ansatz state to compute the energy density as well the thermodynamic potential for this 
system. We minimize the thermodynamic potential to determine the the ansatz functions and the resulting mass gap 
equations. We discuss the results of the present investigation in section IV. Finally we summarize and conclude in 
section V. For the sake of completeness we derive the spinors in the presence of a magnetic field and some of their 
properties, which are presented in the appendix. 



II. THE ANSATZ FOR THE GROUND STATE 

To make the notations clear, we first write down the field operator expansion in the momentum space in the presence 
of a constant magnetic field B in the z— direction for a quark with a current mass m and electric charge q. We choose 
the gauge such that the electromagnetic vector potential is given as A^x) = (0, 0, Bx, 0). The Dirac field operator 
for a particle is given as [2f| 
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^( x ) = 2^ / dl \ )?7 r (»,I\,n) +?r(n,-Pi t )K(»,-R,,n)] e'VV 



(1) 



The sum over n in the above expansion runs from to infinity. In the above, = (p y ,p z ), and, r = ±1 denotes 
the up and down spins. We have suppressed the color and flavor indices of the quark field operators. The quark 
annihilation and antiquark creation operators, q r and q r , respectively, satisfy the quantum algebra 



{«r(n,R,),^«p()} = {qr{n,p^),ql,{n',p'J} = <W'<W<S(l\ - p£)- 



(2) 



In the above, t/ r and 14- are the four component spinors for the quarks and antiquarks respectively. The explicit 
forms of the spinors for the fermions with mass m and electric charge q are given by 



U t (x,p ,n) 
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Vf(x,-n,n) 
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V^MS (9(g)7 n _! - 9(-g)J n _!) 
-i(e„ + to) (9(g)7„ + 9(-g)I n _i) 




(3a) 



(3b) 



(3c) 



(3d) 



In the above, the energy of the n-th Landau level is given as e„ = y/mf + pi + 2n\qi\B = \Jm\ + |p?|. In Eq.s ((3|), 
the functions I n 's (with n > 0) are functions of £ = |gjS|(a; — p y /\qiB\) and are given as 



(4) 



where, -ff n (£) is the Hermite polynomial of the nth order and i_i = 0. The normalization constant c„ is given by 



The functions I n (Q satisfy the orthonormality condition 



(5) 



so that the spinors are properly normalized. The detailed derivation of these spinors and some of their properties are 
presented in the appendix. 

With the field operators now defined in terms of the annihilation and the creation operators in presence of a constant 
magnetic field, we now write down an ansatz for the ground state taken as a squeezed coherent state involving quark 
and antiquarks pairs as [3, [H, HtJ 



\n) =Uq\o). 



(6) 
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Here, Uq is an unitary operator which creates quark-antiquark pairs from the vacuum |0). Explicitly, the operator, 
Uq is given as 



W Q =expf^ J ^g; t (n,2^)< )S (n,^)/ i (n,^)g(n,-^)-/i. 1 



(7) 



where, we have explicitly retained the flavor index i for the quark field operators. In the above ansatz for the ground 
state, f l (n,p z ) is a real function describing the quark antiquark condensates related to the vacuum realignment for 
chiral symmetry breaking. In the above equation, the spin dependent structure a % r s is given by 



1 

I Pi 



-y/2n\qi\B5 r , s - ip z S r ,- 



(8) 



with |pi| = \Zp| + 2n\qi\B denoting the magnitude of the three momentum of the quark/antiquark of i-th flavor (with 
electric charge qi) in presence of a magnetic field. It is easy to show that, aa^ = I, where I is an identity matrix in two 
dimensions. The ansatz functions fi(n,p z ) are determined from the minimization of thermodynamic potential. This 
particular ansatz of Eq.© is a direct generalization of the ansatz considered earlier [23|, [25[, to include the effects of 
magnetic field. Clearly, a nontrivial fi(n,p z ) breaks the chiral symmetry. Summation over three colors is understood 
in the exponent of Uq in Eq. © . 

It is easy to show that the transformation of the ground state as in Eq.© is a Bogoliubov transformation. With 
the ground state transforming as Eq.©, any operator O in the |0) basis transforms as 



o' = u Q ou Q 



(9) 



and, in particular, the creation and the annihilation operators of Eq.© transform as 



€(«)-p*) 



= U c 



Qs{n, — Pj,) 
cos |/| -Or,* sin |/| 
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air sin |/| cos |/| 



9r(n,pJ 



(10) 



which is a Bogoliubov transformation with the transformed 'primed' operators satisfying the same anti-commutation 
relations as the 'unprimed' ones as in Eq.©. Using the transformation Eq. we can expand the quark field 

operator ^( x ) in terms of the primed operators given as, 



(11) 



with q'\Q) = = <ft|0). In the above, we have suppressed the flavor and color indices. It is easy to see that the 
'primed' spinors are give as 



U' r {x,n,p^) = cos\f\U r {x,n,j{.) - al tS sm\f\V s (x,n, -p^) 
Vr{x,n,-p k ) = cos|/|K(z,n, -p^) + a Syr sm\f\U s (x,n,p^). 



(12a) 
(12b) 
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Explicit calculation, e.g. for positive charges yield the following forms of the primed spinors: 



U\(p kl n) = 



y / 2e n (e n + m) 





a 2 p z I n 



VU-PL,n) 



v /2e„(e„ + m) 



v /2e„(e n + m) 



\/2e„(e„ +m) 



-ia,2^2n\q\BI n _i 


ia 2 y / 2n\q\BI„ 
-a2P z I n -i 



ft 2 

ia 2 p z In-i 


iailn-i 

a 2 \/2n\q\BI n - 1 

-iaxl n 




where the functions, a\ and a 2 , are given in terms of the condensate function f(p z ,n) as 

a\ = (e„ + m)cos|/(n,p z )| + |p»| sin \f(n,p z )\ 



a 2 = cos \f(n,p z )\ 



■ shx\ f(n,p z 



(13a) 



(13b) 



(13c) 



(13d) 



(14) 
(15) 



Let us note that with Eq. (|lip . the four component quark field operator gets defined in terms of the vacuum structure 
for chiral symmetry breaking given through Eq. (jSJ) and Eq. (J7| [H, [H[ . 

To include the effects of temperature and density we next write down the state at finite temperature and density 
through a thermal Bog oliubov transformation over the state |f2) using the thermo field dynamics (TFD) 
method as described in Ref.s [3 OL 1 3 1| . This is particularly useful while dealing with operators and expectation values. 
We write the thermal state as 



\QH3,n))=Upjn)=UpM0), (16) 

where Up^ is given as 



with 



^(AM)= / [E / ^ g ;(n,fc 2 )t^(fc z ,n,/3,^(n,^)t + q;(n,fc z )0 + (fc z ,n,/3,/i)£(n,fc^ 



(17) 



In Eq. ()17p . the underlined operators are the operators in the extended Hilbert space associated with thermal doubling 
in TFD method, and, the ansatz functions 0±(n, k z , fl, p) are related to quark and antiquark distributions as can be 
seen through the minimization of the thermodynamic potential. In Eq. (117l) we have suppressed the color and flavor 
indices in the quark and antiquark operators as well as in the functions 9zp. 

All the functions in the ansatz in Eq. (|16|) are obtained by minimizing the thermodynamic potential. We shall carry 
out this minimization in the next section. However, before carrying out the minimization procedure, let us focus our 
attention to the expectation values of some known operators to show that with the above variational ansatz for the 
'ground state' given in Eq.(|16jl these reduce to the already known expressions in the appropriate limits. 

Let us first consider the expectation value of the chiral order parameter. The expectation value for chiral order 
parameter for the i-th flavor is given as 
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(Q03, M )|^iMM) = - £ Nc ^2nV n I — [ m « cos2 /» + IP<I s ' m2 M (1 ~ sin 2 0L - sin 2 0;) 

n=0 ^ ' J €m 

= -r (is) 

where, a n = (2 — <5 n o) is the degeneracy factor of the n-th Landau level (all levels are doubly degenerate except the 
lowest Landau level). As we shall see later, the functions sin 2 8 T will be related to the distribution functions for the 
quarks and antiquarks. Further, for later convenience, it is useful to define cos</> = mi/e ni and sin<% = |p;|/e„i so 
as to rewrite the order parameter Ii as 

h = E f dp z cos (1 - sin 2 tf_ - sin 2 6\) , (19) 

71 = ' n ' 

where we have defined <fi l — </> — 2 As we shall see later, it is convenient to vary <fr l as compared to the original 
condensate function /, given through Eq.([7|). This expression for the order parameter, i^, in the limit of vanishing 
condensates (fi=0) reduces to the expression derived in Ref. [26|. Further, the expression for the chiral condensate 
in the absence of the magnetic field at zero temperature and zero density becomes 

(ij «y ) = -h = - / dp cos $ , (20) 



(2tt) 3 d 

once one realizes that in presence of quantizing magnetic field with discrete Landau levels, one has [3 



f dp \qB\ ^ r 



(2tt)3 (2tt)2 ^ 

This expression for the condensate, Eq.(20) is exactly the same as derived earlier in the absence of the magnetic field 

The other quantity that we wish to investigate is the axial fermion current density that is induced at finite chemical 
potential including the effect of temperature. The expectation value of the axial current density is given by 

(if) = (#7 3 7 5 ^)- 

Using the field operator expansion Eq. flTTl) and Eq. (|13dp for the explicit forms for the spinors, we have for the i-th 
flavor 



(21) 



Integrating over dp y using the orthonormal condition of Eq.([5]), all the terms in the above sum for the Landau levels 
cancel out except for the zeroth Landau level so that, 

Of) = ^# / ** t sin2 9 ~ sin2 «?] ■ (M) 



(27T) 2 

which is identical to that in Ref. [33| once we identify the functions sin 2 6^ as the particle and the antiparticle 
distribution functions for the zero modes (see e.g. Eq. (|33l) in the next section). 

III. EVALUATION OF THERMODYNAMIC POTENTIAL AND GAP EQUATIONS 

As has already been mentioned, we shall consider in the present investigation, the 3-fl.avor Nambu Jona Lasinio 
model including the Kobayashi-Maskawa-t-Hooft (KMT) determinant interaction. The corresponding Hamiltonian 
density is given as 

8 

H = i>\-ia ' V - qBxa 2 + 7 °m)V - G s E [(V^V) 2 - (V^ 5 A A V) 2 ] 

+ K [det f [^{l + 75 )V] + detail - tsM] (23) 
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where i\) % ' a denotes a quark field with color 'a' (a = r 1 g 1 b), and flavor (i = u, d, s), indices. The matrix of current 
quark masses is given by rfi=diag/(m„, m^, m s ) in the flavor space. We shall assume in the present investigation, 
isospin symmetry with m u =m ( i- Strictly speaking, when the electromagnetic effects are taken into account, the current 
quark masses of u and d quarks should not be the same due to the difference in their electrical charges. However, 
because of the smallness of the electromagnetic coupling, we shall ignore this tiny effect and continue with m u = ma 
in the present investigation of chiral symmetry breaking. In Eq. (|23[) . X A 7 A = I, • ■ ■ 8 denote the Gellman matrices 

acting in the flavor space and A = y| If, 1/ as the unit matrix in the flavor space. The four point interaction term 

~ G s is symmetric in SU(3)v x SU(3)a x U(1)v x U(1)a- In contrast, the determinant term ~ K which for the case 
of three flavors generates a six point interaction breaks U(1)a symmetry. In the absence of the magnetic field and 
the mass term, the overall symmetry in the flavor space is SU(3)v x SU(3)a X U{l)v- This spontaneously breaks to 
SU(3)v x U(l)y implying the conservation of the baryon number and the flavor number. The current quark mass 
term introduces additional explicit breaking of chiral symmetry leading to partial conservation of the axial current. 
Due to the presence of magnetic field on the other hand the SU(3)v symmetry in the flavor space reduces to to 
SU(2)v x SU(2)a since the u quark has different electric charge compared to d and s quarks j34|. 
Next, we evaluate the expectation value of the kinetic term in Eg. (1231) which is given as 



T = (n(/3,tx)\ipf(-ia ■ V - qi Bxa 2 WMj3, p)). (24) 



To evaluate this we use Eq. (fTTI) and the results of spatial derivatives on the functions I n (0 (£ = y/\qi\B( 
Py/(h\B))). 

dl r 



x — 



dln-l 

dx 



= v/p [-$/ n _i + y/2(n - l)J„_ a ] ■ (25) 
Using above, a straightforward but tedious manipulation leads to the expression for the kinetic term as 

T = -EE ^CjS^ 1 /* z Kcos0 l + |p I |sin^)(l-sin 2 ^-sin 2 ^) (26) 

The contribution from the quartic interaction term in Eq. (|2"3")l . using Eq. (fT5)l turns out to be, 

V s = -G a <nG8,/*)| £ [(4>\ A 4>? ~ (^7 5 A A ^) 2 ] = -2G S £ p\ (27) 

A=0 i=l,3 

where we have used the properties of the Gellman matrices J2a=o ^tj^M = ^il&jk- 

Finally, the contribution from the six quark interaction term leads to the energy expectation value as 

V det = +K{det f $(l+n)1>]+detS(l-<K)1>]) = -IKhhh- 
The thermodynamic potential is then given by 

n = T + V s + V det -^2fiiPi- (28) 
i=l " 

In the above, [ii is the chemical potential for the quark of flavor i. The total number density of the quarks is given by 

p = E />» = E<*ta> = E E Nc ( 2 n })f l I d p* ^ e - - sin2 +] • ( 29 ) 

1=1 i 71 = i \ ) J 

Finally, for the entropy density for the quarks we have (30l | 

s = _YY N cUn\q l \B r( sin 2^ msm 2^ +cos 2^ lncos 2^) + / (3Q) 



8 



Now the functional minimization of the thermodynamic potential Q with respect to the chiral condensate function 
fi{Pz) leads to 

m. l +AGl l + 2K\e ljk \I J I k M t 

COt0,; = j— : J - J — = — 31 

We have defined, in the above, the constituent quark mass Mj for the i-th flavor as 

M = irii -\- 4GI t + 2K\e l]k \I ] I k (32) 
Finally, the minimization of the thermodynamic potential with respect to the thermal functions 9± (k) gives 

sin2 9 ± n = TIT, ~T ; . (33) 

where, Ui >n = \J Mf + p 2 z + 2n\qi\B) is the excitation energy with the constituent quark mass Mj. 

Substituting the solution for the condensate function of Eq. (f3"Tj) and the thermal function given in Eq. (|3"31 back 
in Eq. (fT9|) yields the chiral condensate as 



ee h = £ Nc \f^ n [ d Pz (Mi) (1 - sin 2 9l_ - sin 2 6\) . (34) 

n=0 



(27T) 



Thus Eq. (f3"2")l and Eq. define the self consistent mass gap equation for the i-th quark flavor. Using the solutions 
for the condensate function as well as the gap equation Eq. (j32|) . the thermodynamic potential given in Eq. (l28t reduces 
to 

N c a n \qi B 

(2tt) 2 

N c a n \qiB 



\ - i\ c u> n \qiJD\ ; 



E 



dp z [ln{l + e-^"*-^} + In {1 + e -/3(^+«)}] 



(2^)2/3 
+ 2G^/ 2 +AKhhh. 

(35) 

The zero temperature and the zero density contribution of the thermodynamic potential (0(T = 0, fi = 0)) in 
the above is ultraviolet divergent, which is also transmitted to the gap equation Eq. (l32[) through the integral Ii in 
Eq. ([34| . In the zero field case (B=0) such integrals are regularized either by a sharp cutoff (a step function in 
|p|) that is common in many effective theories like NJL model (35l-l37| although one can also use a smooth regulator 
|38l KM |49| . The choice of the regulator is a part of the definition of the model with the constraint that the physically 
meaningful results should not eventually be dependent on the regularization prescription. A sharp cutoff in presence 
of the magnetic field suffers from cutoff artifact since the continuous momentum dependence in two spatial dimensions 
are being replaced by a sum over discretized Landau levels. To avoid this, a smooth parametrization was used in 
Ref. [3] in the context of chiral magnetic effects in Polyakov loop NJL model. In the present work however we follow 
the elegant procedure that was followed in Ref. [22| by adding and subtracting a vacuum (zero field) contribution 
to the thermodynamic potential which is also divergent. This manipulation makes the first term of Eq. (l3"51) acquire 
a physically more appealing form by separating the vacuum contribution and the finite field contribution written in 
terms of Riemann-Hurwitz £ functions as 

1 " N c a n \ qi B\ 



EE 



i=l n=0 y ' 

2N, ' 



dp z Jpl + 2n\ qi \B + Mf 



(2 



&E f d P^P 

' i=i J 



j=i 



(36) 
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where, we have defined the dimensionless quantity, Xi = 2 \ q i B\ ' i' e ' ^ ne mass parameter in units of the magnetic field. 
Further, £'(— l,x) = d((z : x)/dz\ z= \ is the derivative of the Riemann-Hurwitz zeta function [40j |. 

Using Eq. ([55|) . the quark antiquark condensate of Eq. ([M|) can also be separated into a zero field (divergent) vacuum 
term, a (finite) field dependent term and a (finite) medium dependent term as 



(ipiipi 



2N C 
(2tt) 3 



dp 



Mi 



VP 2 + M? 



N c Mi\ qi B\ 



E 

n=0 



N c \qi\Ba n 
(2tt)2 



1 Xi 

Xi(l - \nxi) + ]nT(xi) + - In 

I Z7T 



dp z 



M, 



:(sin 2 9 i 



P 4- P 4- P 
vac ' field ' medi 



(37) 



where, we have denoted the three terms in above equation as I vac , I field and I me d respectively. The zero field vacuum 
contributions in Eq. (j36|) as well as in the Eq. (|37|) . can be calculated with a sharp three momentum cut off as is usually 
done in the NJL model (35l. |36|. Thus, e.g., the vacuum part of the order parameter I vac becomes 



2tt 2 



Aa/A 2 +M? - M/loe 



A + 0V 2 + M'l 

Mi 



(38) 



However, since in presence of magnetic field, |p| 2 = p 2 z + 2n\qiB\, the condition of sharp three momentum cutoff 
translates to a finite number of Landau level summations in Eq. (|37[) or in Eg. (1361) with n rnax , the maximum number 

when the component of the momentum in the 



2|g,|B 



of Landau levels that are filled up being given as n max = Int 
z-direction p z = 0. Further, for the medium contribution I me d, this also leads to a cut off for the magnitude of \p z \ 
as A' = ^/A 2 — 2n\qi\B for a given value of n. 

Thus the thermodynamic potential, given by Eq.(|35j). can be rewritten as 



n(/3, /i, b, Mi) = n vac + n field + a med + 2G^] / 2 + ^Khhh 



(39) 



where, 



*2 , 



- 2 ^E / 

i J \v\ 



<A 



dp 

(2^)3 



M? 



Nr 



E 



The field contribution to thermodynamic potential is given by 

3 



field 



2tt 



C'(-l,a;i) - -(a; 2 - aj^lnxj 



(A 2 + M?) X / 2 (2A 2 + Af 2 ) - M* loj 



A + 0Y 2 + Af 2 



a? 
4 



(40) 



(41) 



The derivative of the Riemann-Hurwitz zeta function £(z, x) at z = — 1 is given by [401 ) 

2tan- 1 y + ylog(l + y 2 ) 



1 , 1 2 1 2, 1 , 2 

— xiosx x H — x logaH logx + a: 

2 5 4 2 5 12 5 



C'(-i,x) : 

The medium contribution to the thermodynamic potential is 



exp(27rxy) — 1 



dy. 



V med = ^ ^0p J d Pz [\n{l + e-^-^ )} + ln{l + e-^+«) }] . 



(42) 



(43) 



It may be useful to write down the zero temperature limits of the integrals Q. me d and I m ed- Let us note that at 
zero temperature the particle distribution function sin 2 = 0(/i« — while the antiparticle distribution function 



10 



sin + = 0. The ^-function restricts the magnitude of \p z \ to be less than p\ max = \JP l f ~ 2n\qi\B, where, p^ — 
\/ (if — Mf is the Fermi momentum of the corresponding flavor. Further, this also restricts maximum number of 
Landau levels n max to n\l_ ax = Int[ 2 ^ / | 5 ]. The contribution arising due to the medium to the chiral condensate then 



reduces to 



if 



Similarly, the contribution from the medium to the thermodynamic potential at zero temperature reduces to 



N 

i n—0 



+ 2n| 9i |fl)log 



y/M? + 2n\ qi \B 



(45) 



In the context of neutron star matter, the quark phase that could be present in the interior, consists of the u,d,s 
quarks as well as electrons, in weak equilibrium 

d — > u + e~ + D e - , (46a) 



and, 



(46b) 



(46c) 



leading to the relations between the chemical potentials Hu,Hd,Hs,HE as 

Hs= Hd = Hu+ He- (47) 

The neutrino chemical potentials are taken to be zero as they can diffuse out of the star. So there are two independent 
chemical potentials needed to describe the matter in the neutron star interior which we take to be the quark chemical 
potential Hq (one third of the baryon chemical potential) and the electric charge chemical potential, \x e in terms of 
which the chemical potentials are given by Hs = Hq — |/i e = Hd, Hu = Hq + |/i e and he = —He- In addition, for 
description of the charge neutral matter, there is a further constraint for the chemical potentials through the following 
relation for the particle densities given by 

2 11 

gPn - gPd - gP S ~ PE = 0. (48) 

The quark number densities pi for each flavor are already defined in Eq. (|29l) and the electron number density is given 
by 

9* = E / *• t sin2 9 - sin2 9 +] , (49) 

where, the distribution functions for the electron, sin 2 9 T = l/(exp(w e =p he) + l)j with w e = \Jp\ + 2n|e|-B. 

To calculate the total thermodynamic potential (negative of the pressure) relevant for neutron star one has to add 
the thermodynamic potential il e due to the electrons, to the thermodynamic potential for the quarks as given in 
Eq. (|35|) . The contribution of the electrons is given by 

fie = E§S / ^[ln{l + e -^-^ ) } + ln{l + e-^+^)}] (50) 

The thermodynamic potential (Eq. ((39])), the mass gap equations Eq. (|32|) . Eq. (|36| and the charge neutrality 
condition, Eq. (|48j) are the basis for our numerical calculations for various physical situations that we shall discuss in 
the following section. 
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FIG. 1. Constituent quark masses as functions of magnetic field at zero temperature and zero density. Fig. 1(a) shows the 
constituent quark masses as a function of the magnetic field when the determinant interaction is taken into account. Fig 1(b) 
shows the same when the determinant interaction term is ignored i.e. K = 0. The solid curves refer to constituent masses of 
u- quark, the dotted curve refers to the constituent mass of the d-quarks while the dot-dashed curve refers to the same of the 
strange quark. 



IV. RESULTS AND DISCUSSIONS 



For numerical calculations, we have taken the values of the parameters of the NJL model as follows. The coupling 
constant G s has the dimension of [Mass] -2 while the six fermion coupling K has a dimension [Mass] -5 . To regularize 
the divergent integrals we use a sharp cut-off, A in 3-momentum space. Thus we have five parameters in total, namely 
the current quark masses for the non strange and strange quarks, m q and m s , the two couplings G s , K and the 
three-momentum cutoff A. We have chosen here A = 0.6023 GeV, G S A 2 = 1.835, KA 5 = 12.36, m q = 5.5 MeV and 
m s = 0.1407 GeV as has been used in Ref . [4 1| . After choosing m q = 5.5 MeV, the remaining four parameters are 
fixed by fitting to the pion decay constant and the masses of pion, kaon and 77'. With this set of parameters the mass 
of 77 is underestimated by about six percent and the constituent masses of the light quarks turn out to be Mi = 0.368 
GeV for u-d quarks, while the same for strange quark turns out as M s — 0.549 GeV, at zero temperature and zero 
density. It might be relevant here to comment regarding the choice of the parameters. There have been different sets 
of parameters by other groups also [H, H3, HU for the three flavor NJL model. Although the same principle as above 
is used, e.g., as in Ref [3jJ, the resulting parameter sets are not identical. In particular, the dimensionless coupling 
KA 5 differs by as large as about 30 percent as compared to the value used here. This discrepancy is due to a different 
treatment of the 77' meson. Since NJL model does not confine, and because of the large mass of the rf meson (m r( '=958 
MeV), it lies above the threshold for qq decay with an unphysical imaginary part for the corresponding polarization 
diagram. This is an unavoidable feature of NJL model and leaves an uncertainty which is reflected in the difference 
in the parameter sets by different groups. Within this limitation regarding the parameters of the model, however, we 
proceed with the above parameter set which has already been used in the study of the phase diagram of dense matter 
in Ref. [43| as well as in the context of equation of state for neutron star matter in Ref. |44j |. 

Let us begin the discussion of the results for the case when the charge neutrality condition is not imposed. In this 
case [i e = and all the quark flavors have the same chemical potential [i q . For given values of fi q , T(— and eB, 
we solve the mass gap equation Eg. (f3"2")) self consistently using the expression Eq. (|3"T| for the order parameter. Few 
comments regarding the evaluation of P med of Eq. (|37p may be worth mentioning. In these evaluations, while considering 



2\ qi \B 



as 



zero temperature and nonzero fi q , the Landau levels are filled up upto a maximum value of n, n max = Int 
already mentioned in the previous section. On the other hand, for all finite temperature calculations, the levels are 
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FIG. 2. Constituent quark masses when charge neutrality conditions are not imposed. Fig. 2-a shows the M u at zero baryon 
chemical potential as a function of temperature for different values of the magnetic field. Fig. 2-b shows the same for the 
strange quark mass M s . Both the subplots correspond to nonzero values for the current quark masses given as m u =5.5 MeV 
and m s =140.7 MeV. 



filled upto the maximum Landau level, so that the error in neglecting the higher Landau level is less than 10~ 5 . We 
also observe that for low temperatures, near the cross over transition temperature, there could be multiple solutions 
of the mass gap equation corresponding to multiple extrema of the thermodynamic potential. In such cases, we have 
chosen the solution which has the least value of the thermodynamic potential given in Eq. (j3T))) . We ensure this by 
verifying the positivity of the second derivative of the thermodynamic potential with respect to the corresponding 
masses. 

We show the constituent masses of the three flavors of quarks as modified by magnetic field at zero temperature and 
zero density in FigJT] The magnetic field enhances the order parameters as reflected in the values of the constituent 
masses M u , M& and M s . Because of charge difference, this enhancement is not the same for all the quarks. For the 
couplings G and K as chosen here, the enhancement factors (M(B) — M(B — Q))/M(B = 0) e.g. for eB = 20mJ are 
about 35%, 24%, 12% for u,d and s quarks respectively. We might mention here that the effect of magnetic field on 
chiral symmetry breaking has been considered in NJL model in Ref.[45| without the KMT determinant interaction 
term. For a comparison, we have also plotted in Fig. 1 b, the constituent quark masses without the determinant 
term as a function of the magnetic field. Clearly, the mass splitting between u and d quarks is much larger when the 
determinant interaction is not taken in to account- e.g.(M u (B) — Md(B))/M u (B = 0) = 57% at eB — 20mJ when 
K=0 while the same ratio is about 11% when KA 5 = 12.36. This behavior can be understood as follows. Whereas 
the magnetic field tends to differentiate the constituent quark masses of different flavors, the determinant interaction 
which causes mixing between the constituent quarks of different flavors tends to bring the constituent quark masses 
together. This results in the splitting between the constituent quarks of different flavors becoming smaller when 
determinant interaction is included in presence of magnetic field. Such a behavior is also observed in Ref.fTll] for the 
case of two flavor NJL model. 

We then show the temperature dependence of the constituent quark masses of the u and s quark for zero chemical 
potential for different strengths of the magnetic field, in FigJ5J The phase transition remains a smooth crossover as 
is the case with zero magnetic field. The effect of magnetic field as a catalyser of chiral symmetry breaking is also 
evident. The chiral condensate and hence the constituent quark masses increase in the temperature regime considered 
here when the magnetic field is increased. In these calculations all the Landau levels as appropriate for the given 
magnetic field have been filled up and the lowest Landau level approximation has not been assumed. The qualitative 
aspects of the phase transition remains the same as in case of zero magnetic field. This result is in contrast to the 
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FIG. 3. Constituent quark masses as functions of p q at T = for different strength of magnetic field. Constituent quark masses 
of d,u and s quarks are plotted in Fig a, Fig b and Fig c respectively. 



linear sigma model coupled to quarks [6j where the usual cross over becomes a first order phase transition in the 
presence of strong magnetic field. We might mention here that our results are similar to that of Ref-fllj for the two 
flavor NJL model. 

Next, we discuss the behavior of constituent masses as baryon number density is increased at zero temperature and 
for different strengths of magnetic field. In Fig 131 we show the dependence of the constituent masses M u , Md and M s 
on the quark chemical potential at zero temperature for three different strengths of magnetic field. For zero magnetic 
field, as the quark chemical potential is increased, a first order transition is observed to take place for the value of 
fJ>q = Mc ~ 0.362GeV. At lower values of the quark chemical potential (p q < fi c ), the masses of the quarks stay at 
their vacuum values and the baryon number density remains zero. At p — /i c , the first order transition takes place 
and the light quarks have a drop in their masses from their vacuum values of about 367 MeV to about 52 MeV. The 
baryon number density also jumps from zero to 2.37po, with p — 0.17fm~ 3 being the normal nuclear matter density. 
Because of the six fermion KMT term, this first order transition for the light quarks is also reflected in dropping of 
the strange quark mass from its vacuum value of 549 MeV to about 464 MeV. 

As the magnetic field is increased, the critical chemical potential for this first order transition consistently decreases 
as may be clear from figures Fig. 3a, Fig 3b and Fig 3c. For eB = 10rn^,eB = 15m J, the corresponding values of 
He are 0.327 GeV and 0.323 GeV respectively. For p < p c , the constituent quark masses increases with the magnetic 
field as may be seen in Fig. 3a where we have plotted the d-quark mass as a function of quark chemical potential. 
E.g. for eB = lOm^., the increase in masses of the u and d quarks are about 45 MeV and 30 MeV respectively while 
for strange quarks the corresponding increase in mass is about 21 MeV as compared to zero field case. Since the p c 
decreases with increase in magnetic field, there are windows in the range of chemical potential where it appears e.g. 
in Fig 3a for the u-quark, that the mass decreases with the magnetic field in the range of chemical potential between 
/i = 323 MeV to p = 362 MeV. In this regime however, the chiral transition already has taken place for eB = 15m^.. 
For eB = 10m^, in this regime of chemical potential although the first order transition takes place, the transition is 
weaker compared to the case of eB — 15m^ in the sense that the constituent quark mass is higher compared to the 
case of eB — 15m J. Finally, for the case of zero magnetic field, the transition is still to take place. After the transition 
the ordering in the masses changes depending upon the filling up of the Landau levels in case of nonzero magnetic 
fields. The kinks in the mass variation correspond to filling up of the Landau levels. This decrease of critical chemical 
potential due to the presence of magnetic field has also been observed in dense holographic matter and is termed as 
inverse magnetic catalysis of chiral symmetry breaking [46[. We, however, observe that although the critical chemical 
potential decreases with magnetic field, the corresponding baryonic density increases. This is clearly seen in Fig. 0] 
where we have shown the baryon number density as a function of quark chemical potential for different strengths of 
magnetic field. While for eB = lOmJ, the critical density p c /po = 2.39 is almost similar to the zero field value of the 
same p c / Po — 2.38, the critical density for eB = 15m^ is substantially larger with p c / Po = 3.62. Similar qualitative 
behavior was also observed in Ref . [22| . 

At finite chemical potential, we also observe oscillations of the order parameter with the magnetic field as shown 
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FIG. 4. Baryon density as a function of quark chemical potential for different strengths of magnetic field. 
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FIG. 5. Oscillation of u-quark mass with magnetic field. We have taken fi q = 380 MeV and T=0. 



in Figj5]for u-quark. We have taken the value of the chemical potential as fi q = 380 MeV and taken the temperature, 
T as zero. This phenomenon is similar to the oscillation of the magnetization of a material in presence of external 
magnetic field, known as de Hass van Alphen effect [13]. As observed earlier, this is a consequence of oscillations in 
the density of states at Fermi surface due to the Landau quantization. The oscillatory behavior is seen as long as 

2\q\B < <\Jf*iq — M% and ceases when the first Landau level lies above the Fermi surface [48L l49l|. 

We then discuss the effects of magnetic held on charge neutral dense matter as may be relevant for the matter in 
the interior of the neutron stars. The thermodynamic potential is numerically computed as follows. For given values 
of the quark chemical potential the electric charge chemical potential [Ie and magnetic field eB, the coupled mass 
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FIG. 6. Constituent quark masses for charge neutral matter. Masses of up quarks (a) and strange quarks (b) as functions of 
the quark chemical potential, [i q at T=0 for different strengths of magnetic field 



gap equations given by Eq. ([32|) are solved, using the expression Eq. ([37l) for the order parameter. The values of the 
electric charge chemical potential \ie are varied so that the charge neutrality condition Eq. (PfB1 is satisfied. The 
resulting solutions are then used in Eq.([M| to compute the thermodynamic potential. In doing so, we also check if 
there are multiple solutions to the gap equation and choose the one that has the least value of the thermodynamic 
potential. In Fig. [6j we show the masses of the quarks as functions of chemical potential for charge neutral matter 
for zero temperature. Let us note that at the transition point, the d-quark number density is almost twice that 
of u-quark number density to maintain charge neutrality as mass of s-quark is much too large to contribute to the 
charge density. For this to be realized, the mass of d-quark should be sufficiently smaller as compared the mass of 
u quark to generate the required difference in the number densities. This, in turn, means that [id should be larger 
than ji u unlike the charge neutral case where all the quarks have the same chemical potential. Numerically, it turns 
out that for zero magnetic field this condition is satisfied when /j,d ~ 393 MeV and fi u ~ 318 MeV as compared to 
the common chemical potential \x c — 362 MeV when charge neutrality condition is not imposed. The corresponding 
masses of d and u quarks are about 61MeV and 80MeV respectively compared to the common mass of 52MeV at the 
critical chemical potential when charge neutrality is not imposed. These values of \i u and fid at the transition point 
corresponds to a electron chemical potential of /i e ~ 75 MeV at the transition while the corresponding quark chemical 
potential at the transition is \i q = 368MeV = (i c , which is slightly higher as compared to the value of /i c = 362MeV 
for the case when such neutrality condition is not imposed. At the transition point for the neutral matter, the number 
density of d quarks is almost twice that of u quarks while the electron number density is three orders of magnitude 
lower than either of the quark number densities. As the magnetic field is increased, the constituent masses for the 
three quarks increase for chemical potential smaller than the critical chemical potential. The first sudden drop of the 
strange quark mass (see Fig. 6b) is related to the drop in the light quark masses through the determinant interaction. 
The kink structure in the strange quark mass for higher magnetic fields can be identifiable with the filling of different 
Landau levels. Further, it is observed that a higher magnetic field leads to a smaller value for fi c . Similar to the case 
of non charge neutral matter, however, the critical density becomes higher for increased magnetic field. This magnetic 
catalysis of chiral symmetry breaking is clearly shown in Fig(?l where masses of up and strange quarks are shown as 
functions of baryon density for zero temperature and for different strengths of magnetic field. 

We then study the effect of magnetic field on the equation of state, i.e. pressure as a function of energy for the 
charge neutral matter. This is shown in Fig. [5] for zero temperature. The effect of Landau quantization shows up 
in the kink structure of the equation of state. For smaller magnetic fields, this effect is less visible as the number of 
filled Landau levels are quite large. Further, it may be observed that as the magnetic field is increased, the equation 
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FIG. 7. Constituent quark masses of u and s quark masses as functions of baryon density in units of nuclear matter density po 
for different strengths of magnetic field at T=0. 



of state becomes somewhat stiffer. Since the zero density constituent quark masses increase with magnetic field, the 
vacuum energy density becomes lower compared to the zero field case. Therefore the starting values of pressure in 
presence of field becomes higher compared to the zero field case as seen in Fig(5] For higher densities when chiral 
symmetry is restored, let us first note that the contribution to the thermodynamic potential due to the magnetic field 
as given in Eq. (|41[) increases with magnetic field. Therefore, it means that one has to have a larger chemical potential 
for lower magnetic field as compared to higher field to have the same energy density. So one would naively expect 
the pressure (P — up — e) with lower field to be higher. However, one has to note that chiral symmetry is restored 
at a lower chemical potential for higher magnetic field as may be clear from Fig[6] and hence the number densities 
can become higher leading to a higher pressure. This is what actually happens for larger energy densities for the two 
fields shown in Figj8] Because of the lower critical chemical potential for the case of eB — lOm^, the masses of the 
quarks are smaller and hence the densities becomes higher compared to the case of eB = 5m* for the same energy 
densities leading to a stiffer equation of state. 

Next, we discuss the effects of magnetic field on hot neutral quark matter. Such a condition is relevant for the 
matter in the interior of the proto neutron stars where the temperatures could be of about few tens of MeV. In Fig|5J 
we show the effect of temperature on the masses of the quarks in the magnetized neutral matter. As may be expected, 
the effect of temperature smoothes the behavior of the masses as functions of the quark chemical potential. The 
corresponding equations of state are also shown in Fig[S] 

Let us note that in the equation of state that we have plotted in Fig|S] and FigfTUl the pressure here corresponds to 
the thermodynamic pressure i.e. negative of the thermodynamic potential given in Eq. (|39[) . However, in presence of 
magnetic field, the hydrodynamic pressure can be highly anisotropic [2fj| . [2ll[50| when there is significant magnetization 
of the matter. The pressure in the direction of the field Pm is the thermodynamic pressure P = —CI as defined in 
Eq. (|39[) . On the other hand, the pressure Pj_ in the transverse direction of the applied magnetic field is given by 
P± = P — MB [2l|. Here, M = —dCl/dB is the magnetization of the system. Using the thermodynamic potential 
expression given in Eq. p9p . this can be written as 



M = M med + M field + M c 

where, M me d , the contribution the magnetization from the medium which at zero temperature is given by 

3Cl me d 



M, 



med 



dB 



« + 2n|g i |S)log 



\H±Vz 



(51) 



(52) 
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FIG. 8. Equations of state for charge neutral matter at zero temperature for different strengths of the magnetic field. 
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FIG. 9. Constituent masses of u and s quarks as functions of \x q for charge neutral matter for different strengths of magnetic 
field at T=40 MeV 



where, we have abbreviated A n = \J M 2 + 2n\qi\B. Mfidd is the contribution from the field part of the thermodynamic 
potential flfi e id given as 
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field 



field 



dB 
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12 24 



(53) 
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FIG. 10. Equation of state for charge neutral matter at T=40 MeV for different magnetic field. 



where, 



2 arctan(g/) + y log(l + y 2 ) 
(exp(2irxy) — 1)(1 — exp(— 2iry) 



(54) 



We might note here that Qfi e id term originates from the effect of magnetic field on the Dirac sea so that we can 
recognize that Mfidd is due the magnetization of the Dirac sea. 

Finally, M c in Eq. (|51l) is the contribution to the magnetization arising from the last two terms of the thermodynamic 
potential in Eq. ([55|) and is given as 



M r 



dB' 



(55) 



Here, |^ is the derivative of the quark condensate (-{ipiij))) with respect to the magnetic field given as 
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where, the contribution from the medium at zero temperature is 
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and, the field contribution from the condensate to magnetization 
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(57) 



where, as defined earlier xi = (Mf /2\qi\B) and ^ (x) = j/m ^ s * ne logarithmic derivative of the Gamma function. 

The resulting magnetization at T=0 is plotted in Fig[TT]for \i u = OAGeV = fid- The magnetization exhibits rapid 
de Hass-van Alphen oscillations. The irregularity in the oscillation is due to the unequal masses of the three quarks 
which are calculated self consistently using the gap equation. Unlike in Ref. 21], the magnetization does not become 
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FIG. 11. Magnetization in units of A 2 as a function of magnetic field. The magnetic field in units of m 2 is plotted in a 
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FIG. 12. The parallel Pm and transverse P± pressures of strange quark matter as functions of the magnetic field in units of 
pressure Po for zero magnetic field. Magnetic field in units of m 2 is plotted in a logarithmic scale. We have taken here T=0 
and fj, u = OAGeV = /id- 



constant even after the all the quarks are in the lowest Landau level. This is due to the effect of contribution of 
magnetization from the Dirac sea which is included here along with the Fermi sea contribution given by M mec ii um . 

The transverse and the longitudinal pressure for the system is plotted in Fig |12l Here, we have taken T=0 and 
fi = 400MeV. The oscillatory behavior of the magnetization is reflected in the transverse pressure. The two pressures 
start to differ significantly for magnetic field strengths of about eB = which corresponds to about 10 18 Gauss. 
Such field induced anisotropy in pressure is qualitatively similar as in Ref. [2l|, where, the anisotropic properties of 
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transport coefficients for strange quark matter were considered. While considering neutron star structure, one has 
to also include the free field energy ^B 2 to the total energy and pressure. This term adds to the parallel and the 
transverse pressure with different signs (20j . This can make the pressure in the transverse direction negative leading 
to mechanical instability [2(|. While studying structural properties of compact astrophysical objects endowed with 
magnetic fields such splitting of the pressure in to the parallel and perpendicular direction need to be taken into 
account as this can effect the structure and geometry of the star. 

Finally we end this section with a comment regarding the axial fermion current density induced at finite chemical 
potential. From Eq.s([25|) and Eq.([33]t. 



exp( x /p2 + M 2 _ M .) exp(v^fTMf + Mi)J' 

Thus, although the lowest Landau level contributes to the above expectation value, because of its dependence on the 
constituent quark mass parameter M;, the effects of all the higher Landau levels are implicitly there in Eq. (l55|) as the 
constituent masses here are calculated self consistently using Eq. ([32|) and Eq. (|37|) . Further, because of dependence 
on the constituent quark mass the axial quark current density expectation value also depends upon the coupling in a 
nonpcrturbative manner [a, [5l| . 



V. SUMMARY 



We have analyzed here the ground state structure for chiral symmetry breaking in presence of strong magnetic field. 
The methodology uses an explicit variational construct for the ground state in terms of quark-antiquark pairing. A 
nice feature of the approach is that the four component quark field operator in presence of magnetic field could 
get expressed in terms of the ansatz functions that occurs for the description of the ground state. Apart from the 
methodology being new, we also have new results. Namely, the present investigations have been done in a three flavor 
NJL model along with a flavor mixing six quark determinant interaction at finite temperature and density and fields 
within the same framework. In that sense it generalizes the two flavor NJL model considered in Ref.[ll| for both 
finite temperature and density. The gap functions and the thermal distribution functions could be determined self 
consistently for given values of the temperature, the quark chemical potential and the strength of magnetic field. At 
zero baryon density and high temperature, the qualitative feature of chiral transition remains a cross over transition 
even for magnetic field strength eB=10mJ . The magnetic catalysis of chiral symmetry breaking is also observed. 

At finite densities, the effects of Landau quantization is more dramatic. The order parameter shows oscillation 
similar to the de Hass van Alphen effect for magnetization in metals. However, in the present case of dense quark 
matter, the mass of the quark itself is dependant on the strength of magnetic fields which leads to a non periodic 
oscillation of the order parameter. Although the critical chemical potential, fj, c , for chiral transition consistently 
decreases with increase in the strength of the magnetic field, the corresponding density increases with the magnetic 
field strength. Imposition of electrical charge neutrality condition for the quark matter increases the value for /i c . 
Since the mass of the strange quark plays an important role in maintaining the charge neutrality condition, this in 
turn affects the chiral restoration transition in quark matter. The presence of nonzero magnetic field appears to 
make the equation of state stiffcr. Further, the pressure could be anisotropic if the magnetization of the matter is 
significant. Within the model, this anisotropy starts to become relevant for field strengths around 10 18 Gauss. While 
considering the structural properties of astro physical compact objects having magnetic fields this anisotropy in the 
equation of state should be taken into account as it can affect the geometry and structure of the star. 

We have considered here quark-antiquark pairing in our ansatz for the ground state which is homogeneous with zero 
total momentum as in Eq.([7]). However, it is possible that the condensate could be spatially non-homogeneous with 
a net total momentum [52| - l54j . Further, one could include the effect of deconfinement transition by generalizing the 
present model to Polyakov loop NJL models for three flavors to investigate the inter relationship of deconfinement and 
the chiral transition in presence of strong fields for the three flavor case considered here[55j. This will be particularly 
important for finite temperature calculations. At finite density and small temperatures, the ansatz can be generalized 
to include the diquark condensates in presence of magnetic field [HI, [H, H3 ■ Some of these calculations are in progress 
and will be reported elsewhere. 
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Appendix A: Spinors in a constant magnetic field 

Here we derive the solutions for the spinors for a relativistic charged particle in presence of an external constant 
magnetic field and the field operator expansion for the corresponding fermion field. We shall take the direction of 
the magnetic field to be in the z-direction. We choose the corresponding gauge held = (0, 0,Bx, 0). The Dirac 
equation in presence of the uniform magnetic field is then written as 

= (a ■ n + /3m)V, (Al) 

where, II = p — qA is the kinetic momentum of the particle with electric charge q in the presence of the magnetic 
field. 

Let us first derive the positive energy solutions U(x) of Eq. ljAlj) . For stationary solution of energy E we choose 
U(x) as 



U(x,t) = 



x( x ) 



e' lE \ (A2) 



where </>(x) and x( x ) are the two component spinors. Substitution this ansatz in Eq. (|Al|) leads to 

x{x) = wrk^' (A3) 

so that eliminating x m favor of cf>, leads to an equation for the latter as 

(E 2 -m 2 )(j){x) = {a -n) 2 (j){x). (A4) 

Noting that (cr • II) 2 = II 2 — qa ■ B. With B = (0, 0, B) and with our choice of the gauge — (0, 0, Bx, 0), the above 
equation reduces to 

(E 2 - m 2 )0(x) = [-V 2 + {qBxf - qB(a 3 + 2x Py )}<j>(x.) (A5) 

Next, recognizing the fact that in the RHS, the coordinates y and z do not occur explicitly except for in the derivatives, 
we might assume the solution to be of the form 

0(x) = e l ^ y+p ' z ^f(x) (A6) 

where, f(x) — f a u a , a = ±1 for spin up and spin down respectively with U\ = ( q J ; it_i = ( ^ J ; so that 173/ = af . 
Using Eq. (|A6[) in Eq. (|A5P we have, 

d 2 



t 2 _i_ 



MO = 0, (A7) 



where, we have introduced the dimensionless variables £ = \J\q\B (x — and a a = — — — \ q \ P B +qBa ■ Eq. (|A7p is 

a special form of Hermite differential equation, whose solutions exist for a a = In + 1, n = 0, 1, 2, .... This gives the 
energy levels as 

E 2 na = m 2 + P 2 z + (2n+l)\q\B-qBa. (A8) 
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The solution of Eg.(|A7]l is 

f a (0 = c n e-^H n {£) = J n (0, 
Where H n (^) is the Hermite polynomial of the nth order, with the normalization constant c n given by 



(A9) 



V n!2 n 0T 

The functions /„(£)'s satisfy the completeness relation 

Y,InWn{e) = \q\B5{H-t,'). 

n 

Further using orthonormality condition for the Hermite polynomials, J n 's are normalized as 



da n (Ol m (0 = V\q\B8 n , m , 



I n (£)'s are seen to satisfy the following relations 

d 



dx 



I n (0 - Vk\B[-ZIn(0 + V2nl n -i(0] 



(A10) 

(All) 

(A12) 
(A13) 



Thus with the upper component 0(x) known from Eq. (IA6l) and Eq. (|A9[) . the lower component x( x ) can be evaluated 
from Eq. (|A3|) using the relations Eq. (|A13[) . This leads to the explicit solutions for the positive energy spinors as 
t/(x, t) = U(n, , x) exp(ip^ • x^, - ie n t) with 



U^(x,n,n) 



U x (x,p ,n) 



v / 2e„(e„ + to) 



v /2e„(e n + m) 



/ (e n +m)(e(g)Z n + e(-g)/ n _ 1 ) \ 


P*(e(g )Jn + e(-g)/ n -i) 

V -» V^M^ (Q(q)I n -i + e(-g)I n ) / 

/ o \ 

(e »+m) (6(g) J n -i + 0(-g)J„) 
i^^P (6(g)I„ - 0(-g)J n _i) 

V -p, (6(g)/ n - 6(-g)7„_x) y 



(Al4a) 



(A14b) 



In the above, we have defined e n = \Jp 2 z + to 2 + 2n\q\B and further, we have defined J_x = while for nonnegative 
values of n, 7„(£)'s are given by Eq. (|A9j) . 

In an identical manner, one can obtain the solutions for the antiparticles and the solution can be written as 
V(x,i) = V{x,Tp k ,n)exp{-iYK i . • x^ + ie n t) with 

/ y2^B(6(g)/„-e(-g)J n _ 1 ) \ 
zp z (e(g)/„_i + 9(-g)7„) 


V i(e n + to) (6(g)J n _i + 9(-g)/ n ) j 

iPz (e (g)J„ + 0(-g)J n _i) \ 
v / 2^|gTB(e(g)/„_ 1 - 0(-g)J n _i) 
-i(e n + to) (6(g)J„ + 0(-g)/ n _i) 




V t (x,-p ,n) = ===== 

v /2e„(e„ + to) 



1 



v /2e„(e„ + to) 



(Al5a) 



(A15b) 



y 



The spinors are normalized as 

J dxU r (x,p^,ri) t t/ s (a;,p i ,, to) = 6„, m <5 r , s = / dxK(a;, , n^Kfa, p^ , to) 



(A16) 
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These spinors are used in Eq.([T]) for expansion of the field operators in the momentum space. 
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